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Abstract 

Background: Plants react to iron deficiency stress adopting different kind of adaptive responses. Tomato, a 
Strategy I plant, improves iron uptake through acidification of rhizosphere, reduction of Fe 3+ to Fe 2+ and transport 
of Fe 2+ into the cells. Large-scale transcriptional analyses of roots under iron deficiency are only available for a very 
limited number of plant species with particular emphasis for Arabidopsis thaliana. Regarding tomato, an interesting 
model species for Strategy I plants and an economically important crop, physiological responses to Fe-deficiency 
have been thoroughly described and molecular analyses have provided evidence for genes involved in iron uptake 
mechanisms and their regulation. However, no detailed transcriptome analysis has been described so far. 

Results: A genome-wide transcriptional analysis, performed with a chip that allows to monitor the expression of 
more than 25,000 tomato transcripts, identified 97 differentially expressed transcripts by comparing roots of Fe- 
deficient and Fe-sufficient tomato plants. These transcripts are related to the physiological responses of tomato 
roots to the nutrient stress resulting in an improved iron uptake, including regulatory aspects, translocation, root 
morphological modification and adaptation in primary metabolic pathways, such as glycolysis and TCA cycle. Other 
genes play a role in flavonoid biosynthesis and hormonal metabolism. 

Conclusions: The transcriptional characterization confirmed the presence of the previously described mechanisms 
to adapt to iron starvation in tomato, but also allowed to identify other genes potentially playing a role in this 
process, thus opening new research perspectives to improve the knowledge on the tomato root response to the 
nutrient deficiency. 



Background 

Iron (Fe) deficiency is a yield-limiting factor for a variety 
of field crops all around the world and generally results 
from the interaction of limited soil Fe bioavailability and 
susceptible genotype cultivation [1]. Iron is an important 
microelement for plant life due to its involvement as 
redox-active metal in photosynthesis, mitochondrial 
respiration, nitrogen assimilation, hormone biosynthesis, 
production and scavenging of reactive oxygen species, 
osmoprotection and pathogen defence [2], 

Under aerated conditions at neutral alkaline pH, the 
soluble Fe concentration in soil solution is very low. To 
cope with Fe shortage plants have developed two 
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strategies for its acquisition. The Strategy I (all higher 
plants except grasses) relies on improvement of Fe 
uptake through acidification of soil solution by excretion 
of protons via a plasmalemma P-type ATPase resulting 
in an increased Fe solubility, reduction of Fe + to the 
more soluble Fe 2+ by a Fe m -chelate reductase and plas- 
malemma transport of Fe 2+ by the activity of a Fe trans- 
porter [3]. Some model plants used to study Strategy I 
are dicots such as Arabidopsis thaliana, Solarium lyco- 
persicum and Pisum sativum [4]. 

Plant responses to Fe deficiency have been recently 
analyzed on the basis of large-scale changes not only in 
transcriptome [5-14], but also in proteome [15-20] and 
metabolome [17]. Results of transcriptome analysis are 
influenced by differences in experimental plans, plant 
species and microarray platforms, and thus difficult to 
compare and be generalized. Notwithstanding this draw- 
back, recently, a set of 92 transcripts that robustly 
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reflect the transcriptional response of Arabidopsis to Fe 
deficiency [21], has been described as the "ferrome" by 
Schmidt and Buckhout [21]. The "ferrome" consists of a 
list of transcripts considered to be involved in the basic 
response to iron deficiency. The ferrome is particularly 
enriched in genes related to heavy metal cation trans- 
port and metal homeostasis. Focusing on tomato, a 
plant often used as a model to study Fe deficiency 
(Strategy I) and a crop of economic importance, no 
information is available at genome-wide transcriptional 
level. Two proteomic characterizations of tomato roots 
in response to 1-week of Fe deprivation showed 23 [15] 
and 15 [16] differentially expressed protein spots respec- 
tively. Modifications in proteome suggest changes in 
energy metabolism, sulfur metabolism, response to oxi- 
dative stress and signal transduction. 

In the present work a genome-wide transcriptional 
characterization of tomato roots in response to Fe defi- 
ciency is presented. This approach allowed indentifying 
97 differentially expressed transcripts involved in the 
responses to the nutritional stress. Transcriptional 
changes, mainly related to positive modulation of glyco- 
lysis, TCA and methionine cycle, suggest that tomato 
roots behave similarly to Arabidopsis under Fe defi- 
ciency. Furthermore, flavonoid biosynthesis and root 
morphological changes are revealed as specific tomato 
responses to Fe shortage. 

Results and discussion 

Responses to Fe-deficiency 

Typical responses of Fe-deficiency [22] were observed in 
tomato plants grown for 14 d in the presence of a low 
amount of Fe and thereafter subjected to 7 d of Fe 
deprivation. The chlorophyll content (SPAD index 
value) was reduced in Fe-deficient plants (Table 1). A 
concomitant increase in root Fe ln -chelate reductase 
activity (Table 1) was also observed with values similar 
to those commonly found in roots of Fe-deficient 
tomato plants [23]. Furthermore, Fe-deprived tomato 
plants developed more lateral roots and showed an 
abundant production of root-hairs (Figures 1 and 2). 



Table 1 Leaf SPAD index values and root Fe -chelate 



reductase activity 



Sample 


SPAD 


Fe'"-chelate reductase (mol g" 1 root FW 


index 3 


h 1 )" 


Fe- 


29.5 ± 0.3 


0.37 ± 0.04 


sufficient 






Fe- 


16.8 ± 0.6 


1.41 ± 0.06 


deficient 







a SPAD index value of fully expanded young leaves was determined using a 
SPAD-502 meter (Minolta, Osaka, Japan); mean and SD using data of the three 
biological replicates. 

b Mean and SD of three biological replicates. 




Figure 1 Shoot and root apparatus of tomato plants grown 
under different Fe-supply conditions Leaf detail of Fe-deficient 
(A) and Fe-sufficient (B) plants. Shoot (C) and roots (E) of Fe- 
deficient plants and shoot (D) and roots (E) of Fe-sufficient plants. 



Comparison of root gene expression profiles in Fe- 
deficient and Fe-sufficient tomato plants 

Differences in root gene expression between Fe-suffi- 
cient and Fe-deficient tomato plants were obtained by a 
genome-wide gene expression analysis using a tomato 
chip developed through Combimatrix technology [24]. 
This chip allows monitoring simultaneously the expres- 
sion of more than 25,000 tomato transcripts. Ninety- 




Figure 2 Root apparatus of tomato plants grown under 
different Fe-supply condition Detail of root apparatus of A) Fe- 
deficient and B) Fe-sufficient plants. 
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seven differentially expressed transcripts between Fe- 
deficient and Fe-sufficient tomato roots (75 up-regulated 
and 22 down-regulated) were identified by Linear Mod- 
els for MicroArray (LIMMA) [25] (adjusted p-value < 
0.05; |FC| > 2). This result obtained using a large-scale 
chip reinforce the idea that plant transcriptional 
response to Fe shortage is based on the modulation of a 
relative small set of genes as previously observed for the 
Arabidopsis "ferrome" [21]. 

Manually curated annotation of the 97 differentially 
expressed transcripts was based on results of BlastP ana- 
lysis against UniProt [26] database (Figure 3; Additional 
file 1) using terms of biological process of Gene Ontol- 
ogy (GO) [27]. Sequence grouping in functional cate- 
gories according to the GO terms revealed that the 
most abundant functional category was "metabolic pro- 
cess" both for up-regulated and down-regulated tran- 
scripts (35% and 45% respectively; Figure 3). Other up- 
regulated transcripts belonged to "establishment of loca- 
lization" (12%) and "cell wall organization and biogen- 
esis" (8%), while for the down-regulated transcripts 
"response to stimulus" (9%) was one of the most repre- 
sentative main functional categories (Figure 3). Only up- 
regulated transcripts are present in the "secondary meta- 
bolic process" category (Table 2). Transcripts encoding 
proteins with no sequence homology to known proteins 
were defined as "no hits found" (12% and 5% for up- 
regulated and down-regulated transcripts respectively), 
while a similar percentage of transcripts showed homol- 
ogy to proteins involved in "unknown" biological pro- 
cess (24% and 23% respectively). A selection of 
differentially expressed and discussed in relation to Fe 
deficiency is reported in Table 2. The up-regulation and 



down-regulation of six differentially expressed tran- 
scripts in response to Fe deficiency were confirmed 
through Real-time RT-PCR experiments (Table 3). 
TC ID, description, Real-time RT-PCR relative expres- 
sion value (Fe-deficient vs. Fe-sufficient) and microarray 
fold change value (Fe-deficient vs. Fe-sufficient). Real 
time data were normalized on the EFla gene and were 
performed 4 times on 3 independent experiments. Real- 
time RT-PCR data are expressed as mean ± SD. 

Sixty-one of 97 transcripts are ascribable to adaptive 
responses to Fe deficiency involving Fe homeostasis, 
metabolic process, oxidative stress responses, root mor- 
phological modification, transport processes, hormone 
metabolism and signaling. The others are hardly related 
to specific role showing homology to protein without a 
specific biological process or lacking homology to 
known protein ("no hits found"). 

Fe homestasis 

Our transcriptional analysis confirmed that roots of Fe- 
deficient tomato plants overexpressed genes involved in 
Fe uptake and reduction, including the transcripts 
encoding IRON-REGULATED TRANSPORTER (IRT) 
[28-30] and Fe m -chelate reductase (FRO) [31]. The 
tomato bHLH protein (encoded by LeFER) plays a role 
in Fe-deficiency responses through the expression of 
these two tomato Fe mobilization genes belonging to 
the Fe uptake systems of the Strategy I plants [29,32,33]. 
Our data show the up-regulation of the FER transcript 
(#24) in Fe-deficient roots, which is in agreement with 
the positive modulation of Fe-uptake-related transcripts 
such as LeFROl (#68), LeIRTI (#39) and a transcript 
encoding a NATURAL RESISTANCE-ASSOCIATED 



B 



GO:0050896; 1% 




GO:0050896; 9% 



No hits found; 5% 



GO:0009987; 1% 




0051234; 1% 



GO:0008150 biological process 
GO:0008152 metabolic process 
GO:0009987 cellular process 
GO:0050896 response to stimulus 
GO:0051234 establishment of localization 
GO:0065007 biological regulation 

GO:0071840 cellular component organization and biogenesis 

Figure 3 Functional categories distribution of differentially expressed transcripts. Distribution in main functional categories according to 
the GO "biological process" terms of the 75 up-regulated (A) and 22 down-regulated (B) transcripts in roots of Fe-deficient relative to Fe- 
sufficient plants. Percentage of transcripts is reported for each functional category. 
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Table 2 List of transcripts modulated in response to Fe-deficiency and reported in the Discussion 



# 


ProbelD 3 


Description 1 " 


UniProtlD c 


TC-ID d 


FC e 


p, value, adj f 




Biological process GO:0008150 










7 


ty — i nmi A on ~>~7 c 

\ L 1 yz / z4_ob 5_3 / _b 


SRC2 


B6SND5 


TC1 92724 


3.02 


0.014 


2 


I Cl92763_24l_40_i 


Putative uncharacterized protein 


C6T3H9 


TC1 92763 


2.1 1 


0.030 


3 


1 L 1 9oo 1 y_o(Jl _o4_X2 


Putative uncharacterized protein 


Q401 27 


TC193319 


-3.35 


0.016 


4 


i l i yt>o4 1 _yyo_Jz>_b 


Putative uncharacterized protein 


Q401 27 


TCI 95341 


-3.33 


0.01 1 


5 


I c I y o I z /y_3 /_b 


D-protein 


Q8VWY8 


TC1 96753 


6.1 1 


0.028 


6 


I LI y/(jyb_ojo_J4_Xz 


NtEIG-E80 protein (elicitor inducible gene product) 


Q9FXS6 


TC 197095 


4.69 


0.024 


7 


1 C 1 98323_947_40_S 


Predicted protein 


B9MWQ1 


TC 198323 


2.51 


0.032 


8 


Tnnmci 1/130 3n c 

i l i yyz53_ i43y_3y_b 


Zinc finger protein, putative 


B9SLY6 


TC 199253 


2.82 


0.029 


9 


1 LZUZ30U_4o/_4U_b 


EF-1 alpha-like protein 


049604 


TC202360 


2.15 


0.024 


10 


T/ — m/ir-71 /in c 

I C2045 / 1 _463_40_b 


AiBEl 12, putative 


B9STJ3 


TC204571 


-2.62 


0.023 


1 1 


Tr — mr-in-7 onn ~}C c 

I L2Ub2U/_8yu_3D_b 


Putative uncharacterized protein 


Q40127 


TC205207 


-3.09 


0.048 


12 


I L/U/zUU_8y3_3j_o 


i /it nnn ? o in inn i 
VI 1 _UUU3o/U/UU 1 


E0CVH7 


TC207200 


2.44 


0.042 


13 


T/ — im/m "7on T7 c 

ICzO/4L)/_/3y_3/_b 


Predicted protein 


B9HQW6 


TC207407 


2.36 


0.042 


14 


TC207665_362_36_S 


Putative uncharacterized protein 


A9PCS8 


TC207665 


-2.39 


0.048 


15 


TC2087l2_583_40_S 


Germin-like protein 


Q5DT23 


TC208712 


2.35 


0.046 


16 


TC208745_692_37_S 


Putative uncharacterized protein 


A5C0F7 


TC208745 


4.84 


0.016 


U 


TC20932l_482_35_S 


VITISV_041870 


A5C9V2 


TC209321 


6.94 


0.01 1 


18 


TC209504_302_40_S 


Hydrolase 


Q4PSL3 


TC209504 


2.78 


0.013 


19 


TC211515_728_35_S 


Amino acid binding protein, putative 


B9RBU8 


TC21 1515 


3.66 


0.014 


20 


TC212954_1137_35_S 


Predicted protein 


B9HZ36 


TC2 12954 


2.02 


0.046 


21 


TC213456_100_34_S 


Putative D-protein 


Q6K482 


TC213456 


4.97 


0.035 


22 


TC214599_1004_37_S 


Nodulin-like protein At2g16660/T24l2 1.7 


Q9SLF1 


TC2 14599 


6.87 


0.042 


23 


TC215994_489_36_S 


Putative aminotransferase, class V family protein 


Q1 H8R9 


TC2 15994 


2.08 


0.036 




Biological regulation GO: 0065007 










24 


TC191806_907_40_S 


BHLH transcriptional 


Q5GA67 


TC191806 


6.56 


0.024 


25 


TC191963_10S8_37_S 


Ferritin 


Q308A9 


TC191963 


-4.71 


0.008 


26 


TC194645_664_37_S 


DNA binding protein 


B9SZX2 


TC 194645 


12.34 


0.015 


27 


TC198138_1325_39_S 


Thioredoxin peroxidase 1 


Q7Y240 


TC1 981 38 


7.06 


0.046 


28 


TC203853_451_35_S 


Heterogeneous nuclear ribonucleoprotein 27C 


B9SSS5 


TC203853 


2.90 


0.028 


29 


TC206202_500_40_S 


Thioredoxin II 


B9RLX0 


TC206202 


-2.47 


0.048 


30 


TC215976_330_38_X2 


Ferritin 


Q308A9 


TC2 15976 


-3.64 


0.005 




Cellular component organization or biogenesis GO:0071840 










31 


NP000231_1302_40_S 


Extensin-like protein Dif54 


Q43505 


NP000231 


8.22 


0.01 1 


32 


TCI 91 669_1 238_40_S 


Extensin-like protein Extl 


Q8VWM5 


TC191669 


5.53 


0.003 


33 


TC204863_245_40_S 


Extensin-like protein Extl 


Q8VWM5 


TC204863 


8.98 


0.013 


34 


TC212258_415_40_S 


Extensin-like protein Dif54 


Q43505 


TC2 12258 


9.13 


0.01 1 


35 


TC212487_279_40_S 


Extensin-like protein Extl 


Q8VWM5 


TC2 12487 


6.44 


0.0004 


36 


TC214133_1133_40_S 
Cellular process GO:00099 


Extensin-like protein DiflO 

87 


Q43504 


TC214133 


5.42 


0.024 


37 


TC207486_546_40_S 


Pollen specific actin-depolymerizing factor 2 


Q8H2B6 


TC207486 


4.81 


0.005 




Establishment of localization GO: 0051234 










38 


NP79645 1 _1 558_36_S 


Root-specific metal transporter 


Q84LR1 


NP796451 


12.00 


0.024 


39 


TC191581_1150_36_S 


Iron-regulated transporter 1 


Q9XFB2 


TCI 91 581 


9.42 


0.013 


40 


TC192292_1560_39_S 


Hippocampus abundant transcript 1 protein 


B9SG70 


TCI 92292 


14.79 


0.0004 


41 


TC200857_1001_40_S 


Ammonium transporter 1 member 1 


P58905 


TC200857 


2.05 


0.043 


42 


TC205660_580_35_S 


Metal tolerance protein 


B9GLJ8 


TC205660 


3.77 


0.049 


43 


TC206149_966_36_S 


Aluminum-activated malate transporter 9 


Q9LS46 


TC206149 


2/11 


0.048 


44 


TC208376_922_36_S 


Oligopeptide transporter, putative 


B9SA63 


TC208376 


17.67 


0.004 


45 


TC215768_1 1 1 1_35_X2 


Aquaporin 


Q8W506 


TC2 15768 


-2.03 


0.023 


46 


TC215874_553_40_S 


Sec14 cytosolic factor, putative 


B9S6A7 


TC2 15874 


3.5 


0.009 
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Table 2 List of transcripts modulated in response to Fe-deficiency and reported in the Discussion (Continued) 



47 


Tl Mi^OOH 1111 DO VI 

1 Lzl 6882_1 1 21_38_Xz 


Hippocampus abundant transcript 1 protein 


B9SG70 


Tl ) 1 f O O ~) 

I Lz 16882 


9.29 


0.003 




A An+/nl~\/^ r\rsy/-f)cc f f)'Of~\f~\ 

IVlcluuOllt piuCcSb kjU.UUUi 


3IDZ, LCJ/fJOf/ Ullli/uLIOll (jU.UU ! jy/0 










AO 


TfinAiiz. "\a 1 j 3 c; c 


Phosphoenolpyruvate carboxylase 


D171 <^A 
Y Zl \ j4 


Tr jn/iT) 1 ; 

1 LZU4ZZZ) 


j. .j i 


n n^o 
u.uby 




TV — )i/io7Q 7nn 1£, c 
1 Lz l4y/8_/UU_bo_b 


Phosphoenolpyruvate carboxylase 


D")71 QA 
r Zl 1 j4 


TC — ) 1 ACi7Q 
\ Lz !4y/8 


3 7 A 

b./4 


n mft 
U.Ubo 




IvIclUUUilL piOLeSS kjvJ.UUUi 


JIJZ, V-UIUUOIIL piOCeSS yJvJ.UUUyujO 










50 


Tf 1 CiA QQ A 1 Q C /I 37 C 

1 L 1 y4b84_ 1 8b4_b/_b 


Cysteine-type peptidase, putative 


DnCTvn 
tsyb 1 AU 


Tf 1 H/1 C O/l 

1 L 1 y4b84 


3.61 


0.01 1 


j / 


Tron^^^n 1^73 An ^ 

I L/UjjjU I J / 7_*+U_J 


Vacuolar processing enzyme 1 


ROM 1T0- 

DZIVl 1 1 U 


1 LZUjjjU 


o 1 ft 

Z. I o 


n 07ft 
u.uzo 


jZ 


Tf — )r\Q1 QA incQ 3Q c 
1 LzU8 1 j4_zUjy_3o_b 


Dnfl imiillrn Df~T DIV/IVI 

beta-amyiase ri_i-DiviYi 


L>y4tuy 


Tf — )n<3 1 r^/l 
I LzU8 I D4 


3 n^ 
-b.Uo 


n n/i 7 
U.U4Z 


C3 

JJ 


Tf — )1 CQ7n lAHZ AC\ c 

I Lz I by/U_b4Ub_4U_b 


Protease Do-like 7 


L/8KY zz 


TC — ) 1 QCi7C\ 

I Lz I by/u 


/i m 
4.U I 


n n^i 
U.UbZ 




AAr>tf>l-\r~\li/- r^rrymrc f f)'f~lflf~\ 

iVIciUUUlIC process kjU.UUUi 


3/jz, i-tiiuiuf iiieiuuoiit process kjvj.'juhhzj/ 










54 


Tf 1 mn/i n 1 C/i c /in c 
1 L 1 yzU4y_ 1 o4b_4U_b 


Sulfate adenylyltransferase 


f\A 3 1 Q 3 

U4J 1 O J 


T/* - 1 mn/i n 

i l i yzU4y 


-2.48 


0.029 


cc 

jj 


Tf"1 0303/1 Q7Q 30. Y7 
1 L 1 yjyj4_0Zo_3O_AZ 


Methylthioribose kinase, putative 


RODVQT 
Dyn 1 oZ 


TC~1 O^Q^/I 

i l i ybyj4 


Z. I z 


n n73 

U.Uzb 


JO 


TriQ/j.3P.n 180 3ft ^ 
i v„ i y^t3ou_ i oz_ju_j 


S-odenosylmethionine-dependent methyitfonsferose, putative 


Dy3Z_30 


Tfi n/non 

1 v_ 1 ¥HDO\J 


Z.U3 


n nzm 


s 7 


Tf"1 O^n^l ftQQ 3 C 

i l 1 yjUjz_oyy_3j_j 


CBL-interacting serine/threonine-protein kinase 1 1 


Lzzybz 


Tf 1 o c; n 3 1 
\ l i ybUbz 


") ft") 
-z.oz 


n mn 

U.UbU 


58 


t/* - momn ino 3 /i v~) 

1 l 1 yy 1 uy_/uo_i4_Az 


Phosphofructokinase, putative 




Tf i no 1 nn 

i l i y8 1 uy 


2.00 


0.036 


59 


ty~i oon7~) 1 1 c-i /in c 

1 l 1 yyy/z_ 1 1 t>4_4u_b 


Fructose-bisphosphate aldolase 


LJ^r YAJ 


i l i yyy/z 


2.53 


0.022 


60 


Tf — mi ^cn ^/ic 3^ c 
1 LzU 1 bbu_o4b_bo_b 


Protein phosphatase-2C 


Uoz4oy 


T/ - — >n 1 3 cn 
I LzU I bbU 


2.1 0 


0.028 


6/ 


TV — im cm c/i 7 /in c 
1 LzU 1 oyz_D4/_4U_b 


Xyloglucan endotransglucosylase/hydrolase 14 


DOD1/I C 


Tf — >n i cm 

i LzU i oyz 


7.29 


0.01 3 


0^ 


Tf >nft3^7 3ft 1 3Q 
1 LZUObb /_j0 1 _jo_j 


Catalase isozyme 1 


rbUZ04 


1 LZUubD/ 


7 70 

-z./ y 


n n/ift 

U.U40 


63 


Tf — >i Tf)7Q nci 3/1 c 

i lz i zy/o_yb i_b4_b 


Phosphoenolpyruvate carboxylase 


Q8S91 5 


TC — ) 1 1070 

1 Lz 1 zy/8 


3.62 


0.030 


6-7 


Tf — > 1 /I 0~)£L 1i^Q/1 C 

1 Lz 1 48zo_ I o84_bb_b 


ATP binding protein, putative 


B9RII2 


Tf — ) 1 A O 

1 Lz 1 4ozo 


1 2.39 


0.01 1 


65 


t/ — 11 / oo7 r c\~i a r\ c 

\ Lz 1 48b /_o9 /_4U_b 


Cytokinin oxidase/dehydrogenase 


f~>\ /n ^ yio 

LJVrM8 


Tf 11 /I037 

1 L2 148:5/ 


4.64 


0.01 5 


00 


TC — )1(^C")0 QO£ 30 C 

i lz i obzy_8yo_by_b 


Peroxidase 7 


OOCV33 

L^ybYbb 


Tf~l 1 AOO 

I Lz I obzy 


O ")C 

y.zb 


n m 1 
U.U I I 


67 


TV — 1 1 sc7i C/n /in c 
Lz 1 Dj / Z_b4z_4U_b 


Avr9/Cf-9 rapidly elicited protein 216 


f\o Af\un 
U84UtU 


TC21 6572 


-2.61 


0.042 




IVlcLuuOIIC piOLtbS UU.UUU' 


3 / jz, uxiuuuon reuutuon kju.uujj i m 










f.O 

Oo 


Ti^lOI/ll") 1~)QA 1Q C 

I L 1 y 1 4 1 z_zzo4_Do_3 


Ferric-chelate reductase 


LiotlVILU 


Tf 1 0 1 A 1 1 

1 L I y 1 4 1 z 


I b.z I 


n nno 
U.UU8 


69 


i l i y i oyb_y_j 1 _j 


Superoxide dismutase 


f\~7\/\/ A A 

U / T r\44 


tc" i m ori3 
I L I y I 8yb 


-2.29 


0.030 


70 


i l i y4 i jy_zUoi_jo_3 


Ferric-chelate reductase 


□ nDi I n 

□yniuz 


Tf" 1 n /i 1 3 n 
I L I y4 1 by 


-1 0.76 


0.023 


77 


TriQZL9?7 7ft /il VI 

\ \~ \ y i fzz/_/o_ i t 1 _az 
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Table 2 List of transcripts modulated in response to Fe-deficiency and reported in the Discussion (Continued) 



91 


TC209134_260_40_S 


No hits found 


TC209134 


30.06 


0.007 


92 


TC209988_335_40_S 


No hits found 


TC209988 


6.82 


0.013 


93 


TC211287_216_38_S 


No hits found 


TC2 11287 


7.63 


0.011 


94 


TC212074_241_40_S 


No hits found 


TC2 12074 


6.44 


0.023 


95 


TC212933_306_37_S 


No hits found 


TC2 12933 


9.59 


0.011 


96 


TC214074_254_40_S 


No hits found 


TC2 14074 


2.22 


0.041 


97 


TC215128_252_35_S 


No hits found 


TC215128 


3.86 


0.029 



a ID of TomatoArray2.0 probes 

description of each transcript. Bold discussed; italics undiscussed. 
c UniProtlD [26] of the first hit obtained by BlastP analysis 
d ID of the TC of DFCI Tomato Gene Index (Release 12.0) [92] 
e Fold change value 
^adjusted p-value 



MACROPHAGE PROTEIN1 (LeNRAMPl; #38) [29,32]. 
A positive modulation of FER and Fe mobilization pro- 
teins (IRT1, NRAMP1 and FROl) was not found in two 
proteomics studies performed in tomato roots grown in 
conditions similar to those used in the present work 
[15,16]. Authors justified these results as related to the 
features of proteomic approach, which was not sensitive 
enough to detect FER and not well suited for mem- 
brane-bound proteins. However, functional characteriza- 
tions of fer mutant proved that FER controls the 
expression of the iron-uptake genes [29,32,33]. A tran- 
scriptional behaviour similar to that described in the 
present work was observed for Ambidopsis orthologous 
genes in Fe-starved roots [5,8,11,34]. 

Our data show a strong down-regulation of another 
ferric chelate reductase (#70). Previous results indicated 
that the same transcript (TC194139 of the Release 12.0, 
corresponding to TC124302 of the Release 9.0) specific 
to the Solarium lycopersicum genome, is only slightly 
regulated by Fe and that its function is not essential for 
Fe uptake [34]. Our results showing a negative regula- 
tion of this FRO transcript in response to Fe-deficiency 
in roots, also quantified by Real-time RT-PCR experi- 
ment, confirm that this gene does not play a crucial role 
in deficiency-induced Fe uptake and suggest the involve- 
ment in other biological process. 

Together with the positive modulation of LeFER, 
LeFROl and LeIRTI, we observed for the first time a 



high up-regulation of another bHLH transcript (#26); 
this result suggests that like in Ambidopsis [35] also in 
tomato plants the response to Fe deficiency through 
FER activity may need the interaction with another 
bHLH protein. BlastP analysis against TAIR database 
[36] using protein sequence obtained from the predicted 
coding sequence of the TC194645 showed the highest 
sequence homology with the protein encoded by 
AtbHLH38 (score: 120; Evalue: 9 E-28; identity: 35%; 
positives: 57%) known to interact with FIT, the Arabi- 
dopsis orthologous of tomato FER [35]. 

The down-regulation of two ferritin transcripts (#25 
and 30) is in line with the negative regulation of ferritin 
genes observed in roots of Fe-depleted Arabidopsis 
plants [5,8,11]. It has been suggested that ferritins can 
be involved in Fe homeostasis [37] with a main role of 
plastidial Fe [11]. Arabidopsis nodulin-like genes have 
been recently described to be putatively involved in Fe 
transport and storage under metal cation sufficiency 
[38]. Transcriptional levels of Arabidopsis nodulin-like 
genes were down-regulated at least until 72-h of Fe-defi- 
ciency, while two other nodulin-like genes were not 
modulated in response to different Fe conditions [38]. 
On the other hand in tomato roots we recorded a posi- 
tive modulation of two transcripts encoding protein 
with a nodulin-like domain (#9 and 22). BlastP analysis 
against TAIR database [36] using protein sequence 
obtained from the predicted coding sequence of the 



Table 3 Real-time RT-PCR validation of a set of genes differentially expressed in microarray analysis 



TC ID 


Description 


Real-time RT-PCR [ratio) 


Microarray 


TC208376 


Oligopeptide transporter, putative 


17.87 ± 4.35 


17.67 


NP796451 


Root-specific metal transporter 


15.11 ± 3.97 


12.00 


TCI 91 581 


Iron-regulated transporter 1 


3.07 + 0.59 


9.42 


TC2 16882 


Hippocampus abundant transcript 1 
protein, putative 


8.69 ± 1.53 


9.29 


TC205660 


Metal tolerance protein 


1 .55 ± 0.44 


3.77 


TC194139 


Ferric-chelate reductase, putative 


-12.14 ± 3.15 


-10.76 
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TC214599 (#22) showed the highest sequence homology 
with the protein encoded by At3g43660 (score: 325; 
Evalue: 2 E-29; identity: 54%; positives: 68%) which is 
one of the two Arabidopsis genes not modulated by iron 
[38]. 

Metabolic processes 
Carbohydrate metabolism 

As observed in Arabidopsis root microarray analyses 
[5,8,11] glycolysis-related genes are positively modulated 
in Fe-deficient roots. These transcriptional data fit well 
with the increased activity of glyceraldehyde-3-phos- 
phate dehydrogenase (GADPH), pyruvate kinase (PK), 
and phosphofructokinase (PFK) recorded in response to 
Fe starvation in cucumber roots [39]. Furthermore, 
increased levels of protein related to glycolysis under Fe 
shortage were recently reported in sugar beat [17] and 
Medicago truncatula [20] roots. All together these evi- 
dences are consistent with the idea of a shift from ana- 
bolic to catabolic metabolism. In the present work, a 
fructose-bisphosphate aldolase (FBP aldolase) (#59) and 
a PFK (#58) transcripts were up-regulated under Fe 
deficiency, further confirming changes in primary meta- 
bolism in response to Fe starvation. The positive modu- 
lation of transcript encoding a PFK, an enzyme 
catalysing a protogenic reaction, supports the role of 
glycolysis in different process such as production of 
ATP and H + for H + -ATPase, reducing equivalents for 
ferric chelate reductase and of phosphoenolpyruvate 
(PEP) [40]. Up-regulation of three phosphoenolpyruvate 
carboxylase (PEPC) transcripts (#48, 49 and 63) agrees 
with results of several proteomic and physiological stu- 
dies on the response to Fe deficiency in tomato 
[15,16,41], sugar beat [17] and Medicago truncatula 
roots [20], showing a positive modulation of proteins 
involved both in glycolysis and TCA cycle. PEPC activ- 
ity, through pyruvate consumption, can keep active the 
glycolytic pathway and give a contribution to the con- 
trol of cytosolic pH [40,41]. In Fe deficiency, starch cat- 
abolism was reported to be enhanced both at 
transcriptional [5] and protein level [16]. Our analysis 
revealed a down-regulation of a transcript encoding a 
protein showing homology to a potato chloroplastic fi- 
amylase (#52) involved into starch degradation in plas- 
tids [42]. Since starch catabolism mediated by this 
enzyme occurs in plastids, the negative modulation of 
fi-amylase could be ascribed to other causes than an 
accelerated glycolysis. 

The positive modulation of a cytochrome C oxidase 
(#76) is in line with the transient induction of electron- 
transport-chain genes observed in Arabidopsis [5] and 
consistent with the enhancement of respiration rate 
observed in cumber [39] and sugar beat [43] Fe-deficient 
roots. This behaviour was interpreted as an attempt to 



increase energy production through oxidative phosphor- 
ylation. However, more recently it has been suggested 
that the increased respiration rate in root segments of 
Fe-starved cucumber plants should not be interpreted as 
an increase in mitochondrial activity but rather as the 
result of an increase in the number of less efficient 
mitochondria and of the induction of different 0 2 -con- 
suming reaction [44]. 
Methionine cycle 

Nicotianamine (NA) is considered to be a key molecule 
for long-distance transport of Fe in plants [45]. Proteo- 
mic analysis of Fe-starved tomato roots [16] showed a 
positive modulation of proteins related to metabolism of 
methionine (e.g. methionine synthase), a precursor of 
nicotianamine. An up-regulation of a transcript encod- 
ing a methylthioribose kinase (MTK, #55), another 
enzyme of the methionine cycle, was observed in our 
transcriptional analysis. A positive modulation of MTK 
transcripts was also recorded in roots of rice, a Strategy 
II plant species, under Fe deficiency [46,47] , although in 
this case, related to the necessity to increase the produc- 
tion of mugineic acid. 

A down-regulation of a sulfate adenylyltransferase 
(ATPS; ATP sulfurylase) gene (#54) was observed in this 
work. However, the opposite was found in Arabidopsis 
[11]. The connection between the sulphur nutritional 
status and capability to respond to Fe shortage has been 
demonstrated in both Strategy I and Strategy II plants 
[30,48]. Recently an increase in methionine content 
related to phytosiderophore synthesis without significant 
changes in ATPS activity was described in Fe-deficient 
barley roots [48]. Interestingly, in a recent study of Med- 
icago truncatula root proteome [20] two enzymes 
related to biosynthesis of cysteine and methionine were 
negatively affected by Fe deficiency. 
Protein turnover 

Response to Fe starvation induced the accumulation of 
gene transcripts related to protein turn-over, including 
protease (#53) and peptidase (#50 and 51) involved in 
protein degradation, and a gene encoding a heteroge- 
neous nuclear ribonucleoprotein (#28) that can act in 
the pre-mRNA metabolism preceding protein synthesis. 

The activity of these genes can be related to molecular 
events controlling plant responses to abiotic stresses 
such as nutritional deficiencies [49]. A general increase 
in protein synthesis was reported as a response to Fe- 
deficiency in cucumber roots [50]; furthermore, protein 
recycling in response to Fe starvation was suggested by 
analysis of expression profiles of soybean [13] and pro- 
teomic changes in cucumber [18] and Medicago tronca- 
tula [20] roots. Transcriptomic data presented here are 
in line with the idea that, under Fe deficiency, N recy- 
cling reactions take place, possibly related to the neces- 
sity of additional anaplerotic source of C and N [18,20]. 
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Secondary metabolism 

Phenolic compounds are reported as components of 
root exudates in Fe-deficient Strategy I plants. These 
molecules are involved in chelation and/or reduction of 
rhizospheric insoluble Fe [51,52]. Recently, phenolics 
have been proposed to selectively influence rhizospheric 
microorganisms and be involved in the reutilization of 
apoplastic Fe [53]. Our results showed an up-regulation 
of a flavonoid-3-hydroxylase {#81) gene and of three 
genes putatively involved into flavonoid glycosylation 
(#82, 83 and 84). Two out of the last three genes 
showed sequence homology to a Cathamnthus roseus 
flavonoid glucoside 1,6-glucosyltransferase catalysing 
1,6-glucosylation of flavonol and flavone glucosides [54]. 
A positive modulation of genes related to general phe- 
nylpropanoid pathway (e.g. genes encoding PAL and 
4CL) was reported in Fe-deficient Arabidopsis roots 
[11]. Data of the present work underline the up-regula- 
tion of transcripts involved in a more specific branch of 
phenolic pathway supporting the idea that Fe-starved 
roots might operate flavonoid secretion into the rhizo- 
sphere in order to promote Fe acquisition [55]. 

Oxidative stress responses 

Many proteins involved in antioxidative defence 
response contain Fe in heme group or coordinated to 
the thiol group of cysteine. The modulation of tran- 
scripts and proteins related to oxidative stress response 
in roots of Fe-starved plants seems to depend on the 
species and the experimental conditions [5,13,17,19]. 
However, catalase (CAT) and peroxidase (POX) activ- 
ities are known to be depressed under Fe deficiency 
conditions [56] in tomato leaves. Our data showed a 
main down-regulation of transcripts encoding thiore- 
doxin (TRX) (#29) and detoxifying enzymes catalase 
(CAT; #62), superoxide dismutase (#69 and 71), perox- 
idase (POX; #73, 78, 85 and 87). Two other POX tran- 
scripts (#66 and 77) and a thioredoxin peroxidase gene 
(#27) were, conversely, up-regulated. At protein level, a 
decrease of a CAT was reported while some POXs 
showed higher levels in response to Fe deficiency in 
tomato roots [15]; a different response of peroxidase 
isoforms has also been reported in sunflower [57]. 
Taken together, data of the present work suggest a dif- 
ferent role in response to nutrient stress condition 
between the thioredoxin (ABB) and POX isogenes and 
are in agreement with previous results obtained in 
Medicago truncatula [20] and sugar beat [17] Fe-defi- 
cient roots. 

A germin-like transcript (#15) was up-regulated 
under Fe deficiency, similarly to what has been 
observed in the tomato root proteome analysis [15]. 
The positive modulation of a germin protein reported 
in this proteomic study was justified hypothesizing its 



role in producing hydrogen peroxide for apoplastic Fe 
reduction or in other stress response on the basis of 
sequence similarity to a Nicotiana attenuata germin 
protein [15,58]. 

Root morphological adaptation 

Morphological modifications in roots of Fe-deficient 
plants are well documented [59,60]. Root hairs prolifera- 
tion and development of transfer cells were described in 
Fe-starved tomato plants [29,61-63]. Enhanced forma- 
tion of lateral roots and root hairs was also recorded in 
our experiment (Figure 1). Extensin proteins seem to be 
involved in this latter phenomenon; in fact, we observed 
a strong positive modulation of LeExtl (#32, 33 and 35), 
LeDiflO (#36) and LeDif54 (#31 and 34). It was reported 
that these three genes encode extensin, a structural pro- 
tein putatively conferring physical characteristics of the 
cell wall [64], and act during root hair formation in 
tomato due their predominant expression in root hair 
cells [64,65]. 

The overexpression of an endo-l,4-P -glucanase (#80) 
and a xyloglucan endotransglucosylase/hydrolase protein 
(#61) transcripts could be involved into the cell wall 
loosening [66] associated to root morphological adapta- 
tion to Fe deficiency, as previously observed in a tomato 
root proteomic analysis [16]. These changes in root 
morphology are also supported by the up-regulation of 
a transcript showing homology to a tobacco actin-depo- 
lymerizing factor (#37) related to the pollen tube elon- 
gation [67]. Here, we report a negative modulation of a 
transcript (#45) showing homology to the tobacco aqua- 
porin PIP2 [68] suggesting a role of the PIP2 tomato 
protein in Fe-related morphological root changes. It has 
been suggested that PIP aquaporins can play a role not 
only in root water uptake but also in root development 
[68]. Transgenic plants exhibiting RNAi of PIP2 aqua- 
porins showed a significant increase in the length of pri- 
mary roots [68]. 

Transport processes 

Among the positively modulated transcripts belonging 
to "transport" functional category, a stronger modulation 
(more than 17 times) of an oligopeptide transporter 
(OTP) gene (#44) was observed under Fe shortage. 
BlastP analysis against TAIR database [36] using protein 
sequence obtained from the predicted coding sequence 
of the TC208376 (#44) showed the highest sequence 
homology with the protein encoded by AtOPT3 (score: 
541; Evalue: E-154; identity: 83%; positives: 90%). Posi- 
tive modulation of AtOPT2 and AtOPT3 has been 
recorded in Fe-deficient Arabidopsis roots [8,11,21,69]. 
The plant members of OTP family have been described 
to have different functions in transport physiology such 
as long-distance sulphur distribution, nitrogen 
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mobilization, metal homeostasis, and heavy metal 
sequestration. OPTs can transport glutathione, metal- 
chelates and peptides [70]. It was hypothesized that 
some plants OTPs are able to transport Fe-chelates and 
Fe-NA suggesting a role of these proteins in long-dis- 
tance metal transport in planta [71,72]. A similar func- 
tion can be hypothesized for the tomato up-regulated 
OPT transcript, due to the previously described positive 
modulation of genes related to NA synthesis (see 
methionine cycle paragraph). 

Another up-regulated gene involved in transport phe- 
nomena in response to Fe-deficiency is LeAMTl {#41). 
This tomato gene, firstly isolated from a root hair 
cDNA library [73] was root-specifically expressed and 
positively regulated by ammonium (NH 4 + ) in root hairs 
[73,74]. This result suggests the presence of a linkage 
between NH 4 + uptake and Fe shortage. It has been 
demonstrated that NH 4 + -dependent rhizosphere acidifi- 
cation can improve Fe availability in the rhizosphere 
[75]. Interestingly, nitrate acquisition is limited under Fe 
deficiency [20,76]. 

Furthermore, favouring ammonium uptake with 
respect to nitrate could reduce competition for reducing 
equivalent between nitrogen and Fe acquisition. This is 
also in agreement with the hypothesized N-recycling in 
Fe-deficient plant roots [20]. 

Fe-deficient tomato roots strongly overexpressed two 
transcripts {#40 and 47) encoding a protein sharing fea- 
tures of major facilitator superfamily (MSF) with puta- 
tive transport activities. However, on the basis of their 
sequence homology it was not possible to hypothesize 
an involvement in transport of a specific metabolite or 
mineral nutrient. 

A positive modulation of a gene encoding a metal tol- 
erance protein (MTP; #42) was also recorded. A similar 
behaviour was described in Fe-deficient Arabidopsis 
roots [8,11] and interpreted in the light of low specificity 
of IRT1 transporters, which can transport different 
metals into the Fe-deficient plants. MTP genes might 
therefore play a role in the detoxification of zinc ions 
taken up absorbed under Fe deficiency conditions [11]. 

Fe deficiency in tomato induced the expression of a 
transcript {#43) encoding a protein showing sequence 
homology to Arabidopsis vacuolar aluminium- activated 
malate transporter (ALMT) 9 [77]. An increase in 
organic acids concentration mainly citric and malic 
acids in plant roots under Fe starvation has been 
reported for many plant species [78]. However, a 
decrease of malate levels was observed in tomato root 
tips after 15-d of Fe deficiency while higher contents 
were recorded in leaves and xylem sap [41]. The tomato 
malate transporter gene might be involved in malate 
fluxes at intracellular level and/or in long-distance trans- 
port in planta. 



Hormone metabolism and signaling 

The role of plant hormones in the regulation of Fe defi- 
ciency responses has been extensively studied [79,80]. 
As for Arabidopsis Col-0 accession, we identified the 
positive modulation of a methionine cycle gene, MTK 
gene, that could be related not only to NA synthesis 
(see above) but also to the recycling of methylthioade- 
nosine during ethylene production [11]. This is in agree- 
ment with the hypothesized involvement of ethylene 
and/or auxin in control of hair root production under 
Fe deficiency [79]. Furthermore, data of the present 
work suggest the involvement of other hormones (such 
as gibberellin and cytokinin) in Fe deficiency response 
in tomato roots. Concerning gibberellin, we recorded 
the up-regulation of two gibberellin oxidase 20 
(GA20OX) transcripts {#72 and 74). They can regulate 
root morphological changes through the synthesis of 
active GA1. Indeed, it has been reported that GAs are 
related to tomato root growth [81]. In addition, it was 
shown that the expression of a tomato SlGA20oxl:GUS 
construct in Arabidopsis localized also in columella of 
secondary roots [82], thus suggesting a role in secondary 
root formation. The same authors also reported that the 
expression of the construct was positively affected in 
cotyledons, hypocotyls and roots by benziladenine, a 
synthetic cytokinin [82]. However cytokinins were 
described as negative regulators of root Fe uptake 
mechanism in Arabidopsis through a root-growth 
dependent pathway [83]. Cytokinins could play a similar 
role also in tomato roots. This hypothesis is supported 
by the up-regulation of a cytokinin oxidase/dehydrogen- 
ase (CKX) gene (#65) under Fe starvation. CKX genes 
are involved in regulation of plant growth and develop- 
ment through the control of the cytokinin concentration 
[84]. 

Focusing on signal transduction, we observed the up- 
regulation of a transcript (#60) encoding a protein with 
sequence homology to the MCP5, a protein phospha- 
tase-2C (PP2C) of Mesembryanthemum crystallinum 
[85]. McMCP5 is expressed in roots and is induced in 
response to salt and drought stresses [85]. The tomato 
PP2C gene could play a role in signal transduction of 
stress nutrient condition as Fe starvation. The positive 
modulation of a gene encoding a SRC2 protein contain- 
ing a C2 domain {#1) suggests a putative involvement of 
Ca as secondary messenger. Focusing on Ca, in our 
experiment we also identified a down-regulated gene 
(#57) in Fe-starved roots showing homology to an Ara- 
bidopsis gene {AtSOS2) encoding a CBL-interacting ser- 
ine/threonine-protein kinase 11. Arabidopsis SOS2 
interacts with the Ca binding protein SOS3 (SALT 
OVERLY SENSITIVE 3), thus controlling K and Na 
homeostasis and the response to salt stress [86,87]. In 
addition a negative modulation of a gene (#67) showing 
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homology to a tobacco Avr/Cf-9 rapidly elicited (ACRE) 
transcript encoding a PP2C was observed [88]. Our data 
confirm the results obtained with proteome analysis of 
tomato roots in response to Fe-deficiency, where 
changes in the levels of proteins involved in signal trans- 
duction were reported [16]. The observed transcriptional 
changes in tomato roots can be the result of the percep- 
tion of nutrient stress condition and the following signal 
transduction. 

Taken together, these results suggest that Fe defi- 
ciency responses are, at least in part, dependent on hor- 
monal balance modifications possibly resulting from 
signal perception and transduction. 

Conclusion 

Ninety-seven differentially expressed transcripts were 
identified comparing root transcriptional profiles of Fe- 
deficient and Fe-sufficient tomato plants. Tomato roots 
respond to Fe deficiency by modulating the expression 
of a number of transcripts similar to the model plant 
Ambidopsis. The comparison of tomato Fe-responsive 
transcript set with the Arabidopsis "ferrome" [21], 
encompassing 92 transcripts that robustly represent the 
response to Fe shortage, confirms the involvement of 
the well know homologous key regulatory elements (e.g. 
bHLHs) controlling the expression of transcripts related 
to Fe uptake and translocation (e.g. IRT and FRO). As 
showed by Arabidopsis "ferrome" [21], tomato roots 
modulate transcripts involved in homeostasis of Fe and 
heavy metal cations (e.g. IRT, NRAMP, MTP, ferritin) 
and others cation (e.g. AMT). Both plant species require 
the up-regulation of transcripts related to glycolysis (e.g. 
PFK) and methionine cycle (e.g. MTK), the latter path- 
way being putatively linked to NA biosynthesis in 
response to Fe deficiency. Fe-NA complexes could be 
transported both in tomato and Arabidopsis plants 
through OPTs during the response to Fe shortage. Here 
we describe, for the first time, the modulation of a spe- 
cific branch of phenolic (flavonoids) biosynthesis in 
response to Fe deficiency. In addition, tomato roots 
seem to be more characterized by root morphological 
adaptation, mainly linked to hair root production, as 
suggested by the strong up-regulation of extensin 
transcripts. 

Therefore, this transcriptional study, while confirming 
evidence coming from proteomic studies, allowed identi- 
fying new putative targets for further functional investi- 
gations on the response to Fe deficiency in tomato 
roots. 

Methods 

Plant material, growth conditions and RNA extractions 

Tomato seedling {Solarium lycopersicum L., cv. 'Mar- 
mande superprecoce' from DOTTO Spa, Italy), 



germinated for 6 days on filter paper moistened with 1 
mM CaS0 4 , were grown for 14 days in a continuously 
aerated nutrient solution (pH adjusted at 6.0 with 1 M 
KOH) as reported by Tomasi at al. [22] with 5 uM Fe 
(Fe-EDTA); thereafter, most of the plants were trans- 
ferred for a further week to a Fe-free nutrient solution 
(Fe-deficient) and some tomato plants were transferred 
for a week to a nutrient solution containing 100 uM Fe- 
EDTA (Fe-sufficient plants) as control. Nutrient solu- 
tions were renewed every three days. The controlled cli- 
matic conditions were the following: day/night 
photoperiod, 16/8 h; light intensity, 220 uE irfV ; tem- 
perature (day/night) 25/20°C; RH 70 to 80%. 

At the end of the growing period (27 days), Fe-defi- 
cient tomato plants clearly showed visible symptoms of 
Fe deficiency yellowing of the fully expanded apical 
leaves and proliferation of lateral roots and root hairs 
and increase in the diameter of the sub-apical root zone. 
24 hours before harvesting, all nutrient solutions (both 
for Fe-deficient and Fe-sufficient plants) were renewed 
and the pH was adjusted to 7.5 with 10 mM Hepes- 
KOH. The pH of the growing medium was adjusted to 
this value to mimic as close as possible the conditions 
that are occurring in Fe-deficiency-inducing soil. How- 
ever, in order to favour an equilibrate development of 
tomato plants growing in the Fe-free nutrient solution, 
the exposure to the pH of 7.5 was limited to the last 
two days. 

Roots of Fe-deficient and Fe-sufficient tomato plants 
(27 d-old) were harvested five hours after the beginning 
of light phase. The collected roots were immediately fro- 
zen in liquid nitrogen and stored until further proces- 
sing at -80°C. The collection was repeated in three 
independent cultivations and the roots from six plants 
were pooled for each treatment. 

Ferric-chelate reduction 

To determine the root capacity to reduce Fe(III)-EDTA, 
accordingly to Pinton et al. [89] roots of a single intact 
(Fe-sufficient or Fe-deficient) tomato plants were incu- 
bated in the dark at 25°C for 60 min in 50 mL of an 
aerated solution containing CaS0 4 0.5 mM, BPDS 0.5 
mM, Hepes-KOH 10 mM (pH 5.5) and 0.25 mM of Fe 
(III)-EDTA. Thereafter, the absorbance of the solutions 
at 535 nm was measured at intervals of 15 min and the 
amount of Fe(III) reduced calculated by the concentra- 
tion of the Fe(II)-BPDS 3 complex formed, using an 
extinction coefficient of 22.1 mM' cm" . 

Microarray analysis 

Transcriptional analysis was carried out using a Combi- 
matrix [24], produced by the Plant Functional Genomics 
Center, University of Verona [90] . The chip (TomatoAr- 
ray2.0) carries 25,789 nonredundant probes (23,282 
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unique probes and 2,507 probes with more than one 
target) randomly distributed in triplicate across the 
array, each comprising a 35-40-mer oligonucleotide 
designed using the program oligoarray 2.1 [91]. The 
source of sequence information included tentative con- 
sensus sequences (TCs) derived from the DFCI Tomato 
Gene Index [92] Release 12.0 and expressed sequence 
tags. Eight bacterial oligonucleotide sequences provided 
by CombiMatrix, 8 probes designed on 8 Ambion spikes 
and 40 probes based on Bacillus anthracis, Haemophilus 
ducreyi and Alteromonas phage sequences were used as 
negative controls. Complete description of chip is avail- 
able at the Gene Expression Omnibus [93] under the 
series entry (GPL13934). 

Total RNA was isolated using the Spectrum™ Plant 
Total RNA kit (Sigma- Aldrich) and quantified by spec- 
trophotometry using NanoDrop™ 1000 Tern Scientific). 
RNA quality was evaluated using Agilent 2100 Bioanaly- 
zer (Agilent). Total RNA (1 ug) was amplified and 
labelled using the RNA ampULSe kit (Kreatech). After 
checking the quantity and quality of aRNA by spectro- 
photometry using NanoDrop™ 1000 (Thermo Scientific) 
and the quality subsequent labelling, 4 ug of labelled 
aRNA was hybridized to the array according to the 
manufacturer's recommendations [24]. Pre-hybridiza- 
tion, hybridization, washing and imaging were per- 
formed according to the manufacture's protocols. The 
array was scanned with an Axon GenePix® 4400A scan- 
ner (MDS Analytical Technologies). 

Analysis of raw data was performed using the open 
source software of the Bioconductor project [94,95] with 
the statistical R programming language [96,97]. Back- 
ground adjustment, summarization and quantile normali- 
zation were performed using limma package [25]. 
Differentially expressed probes were identified by linear 
models analysis [25] using limma package and applying 
Bayesian correction, adjusted p-value of 0.05 and a |FC| 
> 2. All microarray expression data are available at the 
Gene Expression Omnibus [93] under the series entry 
(GSE31112). Genes were grouped in main functional 
categories according to the "biological" terms of the 
Gene Ontology [27] assigned to each tomato TC or EST 
(Release 12.0) on the basis of the results of BlastP analysis 
[98] against the UniProt database [26] (Additional file 1). 
Genes without significant BlastP results were classified as 
"no hits found" (Evalue < le-8; identity > 40%). 

Real-time RT-PCR experiments 

0.5 ug of total RNA (checked for quality and quantity 
using a spectrophotometer NanoDrop™ 1000 (Thermo 
Scientific), followed by a migration in an agarose gel) of 
each sample was retrotranscribed using 1 pmol of Oligo 
d(T)23VN (New England Biolabs, Beverly, USA) and 10 
U M-MulV RNase H for 1 h at 42°C (Finnzymes, 



Helsinki, Finland) following the application protocol of 
the manufacturers. After RNA digestion with 1 U RNase 
A (USB, Cleveland, USA) for 1 h at 37°C, gene expres- 
sion analyses were performed by adding 0.16 ul of the 
cDNA to the realtime PCR complete mix, FluoCycleTM 
sybr green (20 ul final volume; Euroclone, Pero, Italy), 
in a DNA Engine Opticon Real-Time PCR Detection 
(Biorad, Hercules, USA). Specific primers (Tm = 58°C) 
were designed to generate 80-140 bp PCR products 
(Additional file 2). Three genes were used as housekeep- 
ing to normalized the data: elongation factor 1-alpha 
EFla (X14449; TC203463; forward: 5'- TGGA- 
TATGCTCCAGTGCTTG-3'; reverse: 5'- 

TTCCTTACCTGAACGCCTGT-3'), histone HI 
(AJ224933; TC192148; forward: 5'- CAAAGGC- 
CAAAACTGCTACC-3'; reverse: 5'-AGGCTTTA- 
CAGCTGCTTTCG-3') and ubiquitin Ubi3 (X58253; 
TC196208; forward: 5'-AGCCAAAGAAGATCAAG- 
CACA-3'; reverse: 5'-GCCTCTGAACCTTTCCAGTG- 
3'). Each Real-Time RT-PCR was performed 4 times on 
3 independent experiments; analyses of real-time result 
were performed using Opticon Monitor 2 software 
(Biorad, Hercules, USA) and R [93] with the qpcR pack- 
age [99]. Efficiencies of amplification were calculated 
following the authors' indications [100]: PCR efficiencies 
were 99.15%, 89.16% and 87.25%, for EFla, HI and 
Ubi3 genes, respectively. The efficiencies for TC191581, 
TC192292, TC194139, TC216882, TC205660, TC208376 
and NP796451 were respectively 94.75, 85.03, 98.92, 
91.55, 92.70, 96.61 and 92.43%. The reported Real time 
data were normalized on the EFla gene. Gene expres- 
sion data were illustrated considering the differences in 
the amplification efficiency of PCR and using the gene 
expression levels in roots of Fe-sufficient plants as refer- 
ence; applying the following formula: 



gene exp ression 



X—f e 



(2 x E x ) 



\c, (x+Fe)— c, {x— Fe)] 



(2 x E y )^ y+Fe ^ c '^ Fe ^ 



Where: E x or y is the percentage value of PCR effi- 
ciency for the amplification of the gene x or y, respec- 
tively; C t (x +Fe ) Ct for the control treatment (+Fe) and 
the considered gene (x); Q(x_ Fe ); Ct for the treated roots 
(-Fe) and the considered gene (x); Q(y +Fe ) Ct for the 
control treatment (+Fe) and the housekeeping gene (y); 
Q(y-Fe) Ct for the treated roots (-Fe) and the house- 
keeping gene (y). 

Additional material 



Additional file 1: Functional annotation of 97 differentially 
expressed transcripts. ProbelD. Fold change, adjusted p-value. reference 
Tentative Consensus sequence in DFCI Tomato Gene Index (Release 12.0) 
http://compbio.dfci. harvard.edu/cgi-bin/tgi/gimain. pl?gudb=tomato. 
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description. UniProtID http://www.uniprot.org/. Biological process GO 
term and E-value are reported for each probe. 

Additional file 2: Primer sequences of Real-time RT-PCR experiment. 

TC ID, description and sequences of forward and reverse primers are 
reported. 
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